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STABLE SOLUTIONS OF SYMMETRIC SYSTEMS ON 
RIEMANNIAN MANIFOLDS 

MOSTAFA FAZLY 


Abstract. We examine stable solutions of the following symmetric system on 
a complete, connected, smooth Riemannian manifold M without boundary, 

-AgUi = ■ ■ ,Um) on M, 

when Ag stands for the Laplace-Beltrami operator, Ui : M ^ R and Hi E 
1 < 2 < m. This system is called symmetric if the matrix of 
partial derivatives of all components of i7, that is M.{u) = is 

symmetric. We prove a stability inequality and a Poincare type inequality 
for stable solutions using the Bochner-Weitzenbock formula. Then, we apply 
these inequalities to establish Liouville theorems and flatness of level sets for 
stable solutions of the above symmetric system, under certain assumptions on 
the manifold and on solutions. 
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1. Introduction 


Suppose that M is a complete, connected, smooth, n-dimensional Riemannian 
manifold without boundary, endowed with a smooth Riemannian metric g = {gij}- 
Consider u = for m € C^(M) that satisfies 

(1.1) — AgUi = Hi{u) on M, 

where Ag stands for the Laplace-Beltrami operator and Hi E for i = 

1,..., m. Let / : M —)■ E be a function in C'^(M). Then, the Riemannian gradient 
and the Laplace-Beltrami operator are given by 


( 1 . 2 ) 

and 

(1.3) 


(Vg/). = 


Agf = diVgiVgf) 


{VWW‘3,f) ■ 


In this article, we frequently refer to the Bochner-Weitzenbock formula that is 

(1-4) + VgAgf ■ Vgf + RiCgi^gf, V g f) . 

Here Uf stands for the Hessian of / as the symmetric 2-tensor given in a local 
patch that is 


(1.5) 


Oifh 
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where is the Christoffel symbol that is 

(1-6) = ]^g^^{dighj + djg,h - dhgij). 

Note that from the definition of the Hessian one can see that 

(1-7) |Vg|Vg/|P < |H/|^ 

The equality holds at p £ M fl {V^/ ^ 0} if and only if there exists /Cfe : M —>• R 
for each fc = 1, • • • , n such that 


( 1 . 8 ) 


Vg(Vg/)fc(p) = ^k{p)ygf{p), 


see 


27, 


for details. Her^ we provide the definition of parabolic manifolds. 
We refer interested readers to 0 2^ 3^ 341 for more information. 


Definition 1.1. A manifold M is called parabolic when for every point p G M. 
there exists a precompact neighborhood Mp of p in M such that for an arbitrary 
positive e there exists a function f^ G that feisi) = 1 for all q G Mp and 

/„|Vg/,|Wg<e. 

For the sake of simplicity, we use the notation djHi(u) = and we assume 

that 

(1.9) diHj{u)djH,{u) > Q, 

for 1 < i, j < m. In addition, Br stands for a geodesic ball of radius i? > 0 centred 
at a given point of M and \Bfi\ is its volume with respect to the volume element 
dVg. We now provide the notion of stable solutions. 

Definition 1.2. A solution u = {ui)i of is said to be stable when there exist 
a sequence of functions f = {fi)i where each Q G C^(M) does not change sign and 
a nonnegative constant X such that the following linearized system holds 

n 

(1.10) - Afi = djH^{u)Q + XQ on M, 

for all i = 1, • ■ ■ ,m. In addition, we assume that djHi(u)QQ > 0 for all 1 < i, j < 

TO. 


For the case of M = M” and scalar equations, that is when to = 1, the above 
notion of stability was derived in connections with the De Giorgi’s conjecture and 
it is used in the literature extensively, see 0, i, El EMI [H S, El . The latter 
conjecture 0 , given in 1978, states that bounded monotone solutions of the Allen- 
Cahn equation must be hyperplane, see 


lia,l2J, 

the case of system of equations, that is when to > 2 


and references therein. For 


(1.11) - Aui = ■ ■ ,Um) on 

the notion of stability 
interested readers to 
originated in the Bose-Einstein condensation and nonlinear optics, 

—Au = Hi{u,v) on R", 

—Av = H 2 {u,v) on R", 


is given in |2Cll - l22l| and references therein. We also refer 
LSlii Eli for the following two-component elliptic system. 


( 1 . 12 ) 


when Hi(u,v) = —uv'^ and H 2 {u,v) = —vu^. For this system, nonnegative mono¬ 
tone solutions, which are Ux„(x)vx„(x) < 0 in R", are of interests. Straightforward 
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computations show that monotone solutions satisfy (jl.lOp for C 2 = 

and A = 0. Note also that dyHi = duH 2 < 0 and dvHiduH 2 > 0. For a similar 
notion of stability, we refer interested readers to [i|, for the Allen-Cahn system 
and to |lj, 0, 21, ^ | for systems with general nonlinearities on bounded 

and unbounded domains. 


We now provide the notion of symmetric systems introduced in 2l| when M = 
R". Symmetric systems play a fundamental role throughout this paper when we 
study system (HI) with a general nonlinearity H{u) = {Hi{u))'^i. 

Definition 1.3. We call system symmetric if the matrix of gradient of all 

components of H that is 

(1.13) 

is symmetric. 

Here is how this article is structured. In Section [2l we provide a stability in¬ 
equality for stable solutions of (11.11) . Applying this stability inequality we prove a 
weighted Poincare type inequality for stable solutions of system dni). In Section [21 
we provide applications of inequalities provided in Section [51 To be mathematically 
more precise, we establish various Liouville theorems regarding stable solutions as 
well as a rigidity result concerning level sets of stable solutions. 

2. Inequalities for stable solutions 

We start this section by an inequality for stable solutions of (HD. For the 
case of a scalar equations, this inequality is used in the literature extensively to 
study symmetry properties, regularity theory, Liouville theorems, etc, regarding 
stable solutions. We refer interested readers to 10, E-El E, E, 111 for more 
information. For the case of system of equations and when M = this inequality 
is given in i 0,112. 


5’ 


Lemma 2.1, Let u denote a stable solution of HI). Then 

m p _ m p 

( 2 . 1 ) ^ / Jd,H,{u)d,H,{u)MjdVg < Yi / 

, Ju Ju 

for any (j) = where (fi € (^^((M) for 1 <i <m. 

Proof. The fact that u is a stable solution implies that there exist a sequence 
C = and a nonnegative constant A such that for all z = 1, • • • , m 


on 


( 2 . 2 ) - AgQ = Yi 

i=i 

Consider a sequence of test functions (f = (^i)™ where 4>i G (^^^(M) for 1 < z < to. 

,2 

Multiplying both sides of (12.21) with ^ and integrating, we get 

(2.3) A / + Y f dgH,iu)QfdVg= f 

Jm JM Si JM Si 


i=i 


Applying the fact that 
(2.4) 


- iv,cy € < 

Si 
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we obtain 

(2.5) / 

Jm si JM 

for each i = 1, • • • , m. Note also that here we have applied the divergence theorem 
for the Laplace-Beltrami operator. Combining (12.51) and (12.31) . we end up with 

n „ ,2 p 

(2.6) ^ \^M^dVg. 

ij=l JfA i=i -'M 

For the left-hand side of (12.61) . straightforward calculations show that 



(2.7) 


m p ,2 ^ p ^2 

= Y. d,H,{u)(:,'fdVg + Y d,H,{u)Q^dVg 

Jm C Jm C 

m „ 

+ Y d^Hi{u)(j)1dVg 
i=iJ^ 

^ /* ±2 m p ±2 

= Y 9,H,iu)Q^dVg + Y d.Hg{u)AVg 


m « 

+ Y diHi{u)(j)1dVg 
i=lJ^ 

m c ( Js2 /A2 \ 

= E / + d,Hg{u)Q^ dVg 

,<g Jm \ C 0 J 

m p 

+ Y d^Hi{u)(j)jdVg 

i=ldM 

m p _ m p 

> 2 ^ / \ldjHi{u)diHj{u)4>i(j)jdVg + Y, / diHi{u)(j)fdVg 

Jm '' Jm 

Y f \JdjHi(u)diHj{u)4)i<t)jdVg. 

/l\/T V 


jj=l 


This finishes the proof. 

□ 


Applying the stability inequality (12.11) . given in Lemma l2.ll we prove a weighted 
Poincare inequality as the following theorem. Note that for the case of scalar 
equations a similar ine quality is established in and used in BE 113 , in the 
Euclidean sense, and in [l6|,ll7| on Riemannian manifolds. For the case of system of 
equations, that is when m > 2, this inequality was derived in in the Euclidean 
sense and later used in 12, 13, El and references therein. 
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Theorem 2.1. Assume that u G C^(M) is a stable solution of fop where m,n > 1. 
Then, the following inequality holds for any rj = G (7^(M) 


m „ 

(2.8) V/ {mCg{VgU,,VgUi) + \'H^,\^ -\Vg\VgU£)rffdVg 

i=l 


+ X1 / ( JdjH,{u)diHj{u)\Vguf\\VgUg\r]^iqg - dgHi(u)VgUi ■ VgUj-qf I dVg 

Jm\'^ j 


< 





Proof. Test the stability inequality m on (fi = \VgUi\irii for the sequence of test 
functions ry = and rji G to get 


(2.9) 


I 


^ [ J djH,{u)d^Hg{u)\VgU^\\VgUj 

Jm 


\ijirijdVg 


< 


E 


/ \Vg{\\/gUM\"dVg =: J. 


We rewrite / as 


m p 

(2.10) ^ = E / \^^H,{u)\\VgU,\^^r]fdVg 

i=i « 

+ E / JdjH,{u)d,Hj{u)\VgU^\\VgUj\rj,r]jdVg. 


A simple integration by parts implies 


m p 

(2.11) J = ^ j WgWgUiW^rjl + |VgUiP|Vg?7iP + 2(77i|Vg'Ui|)Vg|VgMi| • VgIJidVg 

m . 

= E/ |Vg|VgIX.||V + |VgW.nVg77.p--7y2Ag|VgW.Pdyg. 

i=l ^ 

We now apply the Bochner-Weitzenbbck formula, given in (El, to have 


m p 

(2.12) J = E/ |Vg|Vgn.||V + |VgU,nVg7y.|2dyg 

i=l 
m p 

“E / + VgAgMj • VgUi77j^ + RiCg(VgUi, VgUi)77j^dVg. 

.^i7m 
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Combining (|2.12|1 . (j2.10ll and (|2.9I1 implies 

m „ 

(2.13) ^ / \d.H,iu)\\VgU.\^r,^dVg 

< / J djHi(u)d^Hj {u)\VgU^ \\VgUj\r]^r]jdVg 

m « 

m « 

“E / + VgAgWj • + RiCg(VgUi, VgUi)7?,^dVg. 

i=l 


On the other hand, for each i = 1,2,...,to differentiating both sides of the i 


th 


equation in CH) and multiplying with VgUiijf we get 


(2.14) - VgUi • VgAgUiTJi = ^ 8 j H i{u)V gU j ■ VgUiiril. 

1=1 

Integrating both sides of the above equation gives 


(2.15) 


i,l=l 

Note that Ji can be rewritten as 


h ■■= “E / '^gUi-VgAgUi'rjfdVg 

m p 

= ^ / djHi{u)VgUj • VgUirjfdVg =: Ji. 


m p m p 

(2.16) Ji=E/ 5*^r*(w)|VgU.|\2dRg + E / dgH,{u)WgUg-VgU,r,^dVg. 

i=l “'M IM 

From p.l6p and p.l5l) we get 

m p m p 

(2.17) E / d^H,{u)\VgU,fr]fdVg = -E / agi7,(u)VgUg • VgU^ry^dRg 

i=l « Jm 

m « 

“E / VgUi • VgAgM^r^^dCg. 
i=l 

Combining (12.131) and (I2.17p . completes the proof. 


□ 


3. Liouville theorems 

In this section, we provide various Liouville theorems for solutions of as 
consequences of the Poincare inequality given in Theorem 12.11 We start with the 
following Liouville theorem. 

Theorem 3.1. Suppose that u = (u^)™]^ is a stable solution of symmetrie system 
foy where the Ricci curvature is nonnegative and Ricg is not identically zero. 
Assume also that one of the following conditions holds 
(i) M is compact. 
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(ii) M is complete and parabolic and |VgUi| G L°°(M) for each i = 1, ■ ■ ■ ,m. 
Then, each m must be constant for i = 1, - ■ ■ ,m. 

Proof. We start with Part (i). If M is compact, then we set 77 ^ = 1 in the stability 
inequality to get 


(3.1) 


E 


/ {RiCg{V gUi,V gUi) + iHuS" 

JM 


|Vg|VgU,|P) dVg 



Note that 


(3.2) \djHi{u)\\VgUi\\VgUj\ > djHi{u)VgUi ■ VgUj. 


This and dSU imply that 


m ^ 

(3.3) V/ [RtCgiVgU,,VgU,) + \-Hu.f -\^g\VgU^f]dVg<0. 

i=l -IM 

Note that the Ricci curvature is nonnegative and the following inequality holds 

(3.4) \Vg\VgU.\\^ <\nuf. 

From this and (IQ) . we get 

(3.5) f RlCg{V gU,,V gUi)dVg = f [ | J ^ _ |Vg|VgU,||2] dVg = 0. 

JM Jm 

This implies that for every p e M fl {S/gUi ^ 0}, we have 

( 3 . 6 ) \Vg\VgU.\\^{p) = \nuA^ip), 
and 

(3.7) RiCg{VgUi,VgUi){p) =0. 

Note that (13.61) implies that 

(3.8) V g{V gUi)k{p) = Kk{p)V gUi{p). 


From the assumptions on RiCg, we conclude that RiCg is positive definite in some 
open subset of M. This and (1X71) imply that for each i, we have \7gUi{p) = 0 for 


28| 


any p in an open subset of M. Therefore, the unique continuation principle 
implies that each Ui must be constant on M. 

We now consider Part (ii). Since M is parabolic, from Definition 11.11 there exist 


Mp and /e such that fj 


M 


|Vg/, 


l^dVg < e. 


Set r]i := f^ in the Poincare inequality 
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L8|) to get 
(3.9) 


III p 

/ {RiCgiV gU^,V gU^) + iHuif - g\V gniW"^) dVg 

i=l "'-^p 

+ i\9jHi{u)\\ygUi\\VgUj\-djHi{u)VgU^-S/gUj)dVg 

JM, 

m p 

< V / {RlCg{VgU,,VgU,)+\nuf -\Vg\VgU,\\^) f,dVg 
i=l "'m 

m ^ 

< E / 

i=l 

m « 

< E|I^9«^II'“(M) / 

i=l -'M 

< einax{||VgMi||^oo(M)}- 

Sending e —5> 0 and using the fact that the inequality (13.21) holds, we get the following 
estimate that is a counterpart of (13.31) on Mp, 

m « 

(3.10) E/ [R^Cg{S/gUi,S/gUi) + \Hui\^ -\^ g\'^gUi\\^]dVg <0. 


i = l "Jp 

Applying similar arguments given for Part (i), completes the proof. 


□ 


The next theorem deals with compact manifolds when Ricg is identically zero. 
We show that when the manifold is compact and RiCg is precisely zero, then stable 
solutions must be constant. Comparing Theorem 13.11 with the following theorem 
may help us better understand the role of RiCg. We refer interested readers to [2^ 
for more information. 

Theorem 3.2. Let u = (Mi)ti ® stable solution for symmetric system 031). 
Suppose that Wl is a compact and connected Riemannian manifold and RiCg is zero. 
Then, each Ui is constant for i = 1, ■ ■ ■ ,m. 

Proof. For each i, differentiate (lEB and multiply with VgUi to get 

m 

(3.11) V gUi • V g^gUi — ^ ) dj (n) V gUj • V gU^ . 

i=i 

Note that dUD implies that for each i 

(3.12) \^g\^gu£ <\nuf. 

Adding (13.111) and (13.121) . we get 

m 

(3.13) + '^gUi ■ VgAgMj > |Vg|VgUi|P - Bj Hi{u)V gUj ' V gUi . 

i=i 
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Since Ricg is assumed to be zero, the Bochner-Weitzenbock formula implies 

(3.14) i Ag| VgMjp = iRui P + ^g^gUi ■ VgU^. 

Combining this and (13.131) we get 

^ m 

(3.15) —Ag|VgUi| > |Vg|VgUi|| ~ ^ ( 9jHi{u)\7gUj ' VgUj. 

i=i 

Integrating this on M and considering the fact that M does not have a boundary, 
for each i, we obtain 


(3.16) 


IVgIVgM.lP ■ VgM* 

i=i 


dVg < 0 . 


Taking summation on index 1 < z < m for (13.161) . yields 

m p 

(3.17) h ■■= y, [\Vg\VgU,W^dVg-d,H,{u)\VgU,\'^]dVg 

i=l 
m p 

< 'y djHi{u)VgUj ■ VgUidVg =: h- 

Jm 

On the other hand, from the stability inequality (EH) for symmetric systems, we 
have 

m „ 

(3.18) hi^) := y [\Vg(l,,\^-d,H,iu)(l,^]dVg 

i=l 
m p 

- y \djHi{u)\(l)i<l)jdVg =: hifi)), 

when (p = ((/ii)™ 1 is a sequence of test functions. Setting pi = |VgMi| in the above 
IsiP) and hiP), we get 


(3.19) 


h=h{P) and l 2 <h{P). 


This implies that p = {pi)’^-^^ when pi = |VgMi| is the minimizer of the energy 
J{p) := h{p) — hip) > 0. Therefore, 

n 

(3.20) - AglVgltil = ^9gi7i('U)|VgMg| OU M, 

i=i 

for each z = 1, • • • , m. From the compactness of M, there exists a point pi G M such 
that Uiipi) = maxMUi. This implies that VgZZi(pi) = 0. From the strong maximum 
principle, we get | VgZZij = 0 for each z = 1, • • • , m. This completes the proof. 

□ 


We now provide the following Liouville theorem in lower dimensions for any non¬ 
linearity H = {Hi)YPi with nonnegative components Hi. Note that for a complete, 
connected, Riemannian manifold M with nonnegative Ricci curvature and dimen¬ 
sion n, the volume of a ball of radius R, denoted |i?R|, is bounded by i?". This fact 
implies that the following Liouville theorem is valid for rz < 4. 
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Theorem 3.3. Suppose that u = (ui)™ i is a bounded stable solution of HJH where 
Hi > 0 for each i = 1, ■ ■ ■ ,m. Let the following decay estimate hold 

(3.21) liminfi?"'‘|Bfl| = 0. 

R—^oo 

Assume also that the Ricci curvature o/M is nonnegative and RiCg is not identically 
zero. Then, each ut must be constant for i = 1, ■ ■ ■ ,m. 

Proof. Suppose that ( G 2, 2]) —>• [0,1] where C = 1 on [—1,1]. For p G M 

and i? > 0, set 

(3.22) Cnip) = C , 

where dg is the geodesic distance. Therefore, Cfi S C))°(M) satisfies Cfl = 1 on Bn 
and Cfi = 0 on M \ B 2 R and 

(3.23) l|VgCfi||L“(B2H\BR) < 

From the assumptions, each Ui is bounded. Multiplying both sides of system (ED 
with (ui — ||Mi||L=o(M))CB from the fact that Hi > 0, we get 

(3.24) ^fi(w)[ui - ||mj||l“(m)]Cb < 0 in M. 

From this and the equation (ED, we have 

(3.25) - - ||mj||l~(m))Cb < 0 in M. 

Doing integration by parts, for each i = 1 • • • , m, we obtain 

(3.26) / \VgUt\^f\dVg <2 [ |VgMi||VCB|(||MilU~(M) - M*)CBdVg. 

J B2R J B2R 

Applying the Cauchy-Schwarz inequality yields 


(3.27) 


[ \VgU.\\ldVg <C f \VgU?dVg 

j B 2 R ^ B 2 R 


where the constant C is independent from R. From the definition of fn and applying 
(13.271) , we get 


m . m . 

(3.28) ^ / \VgU,\^dVg < E / \^gn,\"CldVg 

i=l dBR JB2R 

< c[ iVgCflpdCg 

JBnR 


i?2 ■ 


We now apply the Poincare inequality (12.81) with the test function rji = Cr. In the 
light of (I3(28L the right-hand side of (12.8p becomes 

(3.29) liminf^ / \VgU^\^\Vgr]^\^dVg < liminf-^^/ \VgU^\^dVg 


' B 2 R 

< Climinf = 0. 
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Therefore, 

m „ 

(3.30) V/ {RiCg{VgU^,VgUr) + \nuA'^-\^a\ygU^\\^)dVg 

i=l 

(3.31) +y^j {\d,H,[u)\\VgU,\\VgUg\-dgH,{u)VgU,-VgU,)dVg<0. 

Jm 

This implies that 

m « 

(3.32) V/ {RlCg{VgU,,VgUi) + \'Hu,\^ -\Vg\VgU,\\^)dVg<0. 

The rest of the argument is very similar to the ones provided in the proof of Theorem 

o 

□ 


The ideas and methods applied in the above proof are strongly motivated by 
the ones provided in 


14j . when M = R”, and in 



when the domain is a 

Riemannian manifold. Note that similar idea are used in [22j| for the case of system 
of equations on K”. Lastly, in two dimensions, we have the following rigidity result 
for level sets of solutions. Note that as it is shown in till, for the case of scalar 


equation, the following flatness result does not hold for the Allen-Cahn equation 
in K" with n > 9 endowed with its standard flat metric. We assume that RiCg 
is identically zero that is equivalent to the Gaussian curvature to be zero in two 
dimensions. We refer interested readers to [3 12, [nisi for similar flatness results. 


Theorem 3.4. Suppose that u = (wi)(Ti is a stable solution of E2P and each 
|VgMi| G L°°(M) where M is o Riemannian manifold with dim M = 2. Assume 
that RiCg is identically zero. Then, every connected component of any level set of 
each Ui, on which |VgMi| does not vanish, must be a geodesic. 


Proof. Note that assumptions imply that M is parabolic and it has nonnegative 
Gaussian curvature. Therefore, for each p £ M there exist a set Mp and a function 
/e such that 


(3.33) 


/ \^gM^dVg<e. 

JM 

Set r]i = /g in the Poincare inequality (|2.8p to get 

m p 

(3.34) V] / (RZCg(VgMj, VgMj) + - |Vg|VgUi|P) dVg 

i=l Jm. 

+ i\9jJJii^)\\'^gU^\\'^gUj\-9jH^{u)ygU^-VgUj)dVg 

. / Jm„ 


< emax{||VgM,|||o„(M[)}. 


When e approaches zero, we have 

(3.35) [ RlCg{W gU,,V gUi)dVg = [ [ | J ^ _ | Vg | Vg U, | | ^ ] dRg = Q . 

Jm„ Jm„ 


Therefore, 

(3.36) 


\^9\'^aU.\?{p) = \nu,\\p), 
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for every p S Mp n {VgUi 0}. Since the equality in p.7p holds, there exists 
Kfe : M —)• R for each k = 1, • • • ,n such that 

(3.37) y g{^ gUi)k(j>) = Kk(p)y gU,{p). 

For each i = 1, • • • ,m, consider a connected component F of {ui = C}ri{VgUi ^ 0} 
that is a smooth curve. Let 7 ^ : R —^ M with 

(3.38) l7zP = l- 
It is sufficient to show that 

(3.39) 7 >(t) + T\^{pn\tn^{t) = 0 in R. 

We show that (13.391) holds for an arbitrary value < S R and p= Differentiating 
(13.381) with respect to t we get 

(3.40) 0 = 

We now use normal coordinates at some fixed point p G A4. Suppose that 

(3.41) gxgiP) = Sxp, dkgxgiP) = 0 and r^^(p) = 0. 

Therefore, 

(3.42) 7 ,(t) 7 ,(t) =0. 

Note that for any point on F we have Ui{j{t)) = C. From this we get 

(3.43) 0 = dxUi{j{t))ii^{t), 
and 

(3.44) 0 = dxf,Ui{"f{t))ii^{t)ii^(t) + dxu^{j{t))'yi^(t). 

Combining (13.441) and (13.371) for p = p and k = p, we end up with 

(3.45) 0 = [Kp{p)it'^ii)][dxUi{p)ii^{t)] + dxUi{p)'yi^{F). 

Substituting (I3.43P in p.45p . we get 

(3.46) 0 = dxui{p)'yi^(t). 

From (13.461) and (13.421) we conclude that 7 i'^(t) = 0, since 7 )'^(t) is orthogonal 
to that is tangent to {ui = C} at p, and to dxUi{p) that is orthogonal to 

{ui = C} at p. This and (13.411) imply that (13.391) holds at i. This completes the 
proof. 

□ 

We end this section with the fact that for the most of our main results in this 
section, we assumed that the Ricci curvature is nonnegative. We would like to refer 
interested readers to where hyperbolic spaces are discussed, and references 
therein for rigidity results when the curvature is negative. 
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